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ABSTRACT 

We investigate the possibility of mapping large-scale anti-cyclonic vortices, 
resulting from a global baroclinic instability, as pre-cursors of planet forma- 
tion in proto-planetary disks with the planned Atacama Large Millimeter Array 
(ALMA). On the basis of three-dimensional radiative transfer simulations, im- 
ages of a hydrodynamically calculated disk are derived which provide the basis 
for the simulation of ALMA. We find that ALMA will be able to trace the theo- 
retically predicted large-scale anti-cyclonic vortex and will therefore allow testing 
of existing models of this very early stage of planet formation in circumstellar 
disks. 

Subject headings: hydrodynamics, radiative transfer — techniques: interfero- 
metric — (stars:) circumstellar matter, planetary systems, pre-main sequence 

1. Introduction 

Klahr & Bodenheimer (2002) show that a global baroclinic instability represents a source 
for turbulence leading to angular momentum transport in Keplerian accretion disks with a 
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radial gradient in entropy. Their hydrodynamical simulations show that this baroclinic flow 
is unstable and produces pressure waves, Rossby waves, and vortices in the r — plane of the 
disk. Most interestingly, these hyper-dens anti-cyclonic vortices form out of little background 
noise and become long-lasting features, which have been suggested to lead to the formation 
of planets (e.g. Adams & Watkins 1995). These pre-planetary matter concentrations have 
a surface density up to four times higher then the ambient medium. They furthermore 
concentrate dust in their centres, stressing their importance for the planetary formation 
process (e.g. Barge & Sommeria 1995). 



2. Hydrodynamic and radiative transfer simulations 

Using the code TRAMP (Klahr et al. 1999), we perform 2D hydrodynamic simulations 
of the evolution of the inner region (1-10 AU) of a circumstellar disk with an initial mass 
of K 2 x 1CT 3 M . The comparison with previous models of circumstellar disks around 
T Tauri stars shows (see, e.g., Wood et al. 2002, Cotera et al. 2001) that this corresponds 
to a total disk mass on the order of 10~ 2 M Q . As described by Klahr & Bodenheimer (2002) 
we treat the evolution of the vertically integrated density £ and internal energy ~ ST cent . 
The initial model had a slope in surface density of i?~ 75 and a temperature slope of R~ 10 
which represents a constant pressure scale height of H/r = 0.055. The vertically integrated 
pressure is then given by the ideal gas equation. Sources Q + for internal energy are the 
usual PdV-work as well as the dissipation of possible shocks via a simple von Neumann- 
Richtmyer viscosity (see Stone & Norman 1992). A sink for the internal energy is a black 
body cooling function Q~ = aT e 4 ff mimicking the effect of optically thick radiation in the disk 
and cooling with the effective (surface) temperature T e & = T cent (4/3r) 1 / 4 with the optical 
depth r as a function of temperature and surface density. The simulation uses 128 cells 
in the azimuthal direction and 128 logarithmically distributed cells in the radial direction. 
The initial distribution was disturbed by a 0.1 percent perturbation in density and left to 
evolve freely afterwards under the influence of the global radial entropy gradient. After 10 4 
yrs a huge anti-cyclonic hyper-dens vortex has formed (see Klahr & Bodenheimer 2002). 
Afterwards the 3D data was re-established by assuming a local vertical pressure equilibrium. 

Based on the density and temperature distributions we performed radiative transfer 
simulations using the three-dimensional continuum radiative transfer code MC3D (Wolf 2002; 
see also Wolf et al. 1999, Wolf & Henning 2000). The goal of these simulations was to obtain 
images at frequencies of 345 GHz and 900 GHz. While the first frequency was chosen in order 
to allow comparison with former investigations focused mainly on this frequency, the second 



- 3- 



frequency marks the planned upper limit of the frequency range to be covered by ALMA 2 
and is of special importance for the observation of small spatial structures. The resulting 
images provide the basis for simulations of the continuum observations with ALMA discussed 
in §3. 

Since we consider continuum observations at submillimeter wavelengths, the circum- 
stellar dust is the dominant source of emission. The dust grains are assumed spherical, 
consisting of a mixture of 62.5 % astronomical silicate and 37.5 % graphite (optical data 
from Weingartner & Draine 2001 3 ). For the graphite we adopt the usual "| — |" approxi- 
mation (Draine & Malhotra 1993). The size distribution n(a) of the grains follows a power 
law, n(a) oc a -3 ' 5 , with grain radii in the range 0.005 jum < a < 1 /im. The gas-to-dust mass 
ratio amounts to 100:1 in our simulations. 



3. Simulations of observations with ALMA 

For simulation of the observations with ALMA we use the simulator software published 
by Pety et al. (2001). Since the goal of the observation is to find small-scale structures in 
circumstellar disks, we consider the ALMA configuration with the longest baselines (max. 
baseline: ^13 km) in which the 64 12m-antennas are distributed in rings. In accordance 
with the suggestions given by Pety et al. (2001), the following main types of errors are 
introduced in order to make the simulations realistic: (a) frequency-dependent receiver- and 
system temperatures as given by Guilloteau (2002), (b) random pointing error during the 
observation with a maximum value of of 0.6" in each direction, and (c) 30° phase noise. 
Further error sources, such as amplitude errors and "anomalous" refraction (due to the 
variation of the refractive index of the wet air along the line of sight) are considered as well. 
The observations are simulated for continuous observations centered on the meridian transit 
of the object. The object passes the meridian in zenith where an opacity of 0.15 is assumed. 
The bandwidth amounts to 8 GHz. 

Fig. 1 shows simulated images of the undisturbed disk in comparison to the disk with 
the evolved baroclinic instability at 345 GHz and 900 GHz for object distances of 50 pc and 
140 pc. The hyper-dens vortex is clearly visible even at the lower frequency. However, in 
the case of an assumed distance of 140 pc, the vortex causes a slight asymmetry in the 
reconstructed image only and could be hardly identified. Since several well-studied nearby 
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star-forming regions which can be observed from the prospective ALMA site are located at 
a distance between 140 pc and 200 pc, such as Taurus, Lupus, Ophiuchus, Corona Australis, 
and Chamaeleon, we concentrate on simulations at 900 GHz in the following. Furthermore, 
even more distant star- forming regions (e.g., in Aquila, Vulpecula, Vela, Perseus, Puppis, 
and Orion) are considered by comparison of simulated 900 GHz observations of the evolved 
disk seen under different inclinations in distances up to 300 pc in Fig. 2. While the main 
characteristics of the vortex (pre-planetary matter concentration), such as its spatial extent 
and brightness contrast to the surrounding medium clearly can be extracted from the images 
of objects at a distance of 140 pc, even the spiral pattern in the circumstellar disk (which is 
another, weaker feature of the disk in this evolutionary state) could be resolved in the case 
of a very close young stellar object such as TW Hya. At a distance of 300 pc the vortex is 
visible at low inclinations of the disk only, while it causes a slight asymmetry of the image 
at higher inclinations. 

Finally, we consider possible observing strategies (selected objects vs. surveys) which 
mainly results in the question of required total integration time per object. In Fig. 3 images 
of the evolved disk seen face-on resulting from total integration times of 5 min to 8 h are 
shown. It turns out that exposure times of 5 min are sufficient to detect the vortex and to 
derive its roughly parameters. In order to map the spiral density structures in the case of 
an object at a distance of 50 pc, an integration time of 1-2 h is required. A longer exposure 
time does not significantly increase the information content of the reconstructed image. 
For objects in nearby star-forming regions (distance: 140 pc), a quick "snapshot" survey is 
sufficient since even much longer integration times will not allow derivation of further details 
of the turbulent structure of the disk. 



4. Conclusions 

We found that ALMA will be able to map the anti-cyclonic hyper-dens vortices from 
the global baroclinic instability representing an early stage of planet formation. In nearby 
star-forming regions, the main parameters of the predicted global baroclinic instability in 
proto-planetary disks can be derived within several minutes of integration time. For even 
closer objects (within a distance of about 50 pc) a detailed structure of the spiral density 
pattern will be possible. We want to stress that the largest planned baselines and high 
frequencies are required for these observations. 
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AS: 315 GHz :900 GHz 




Fig. 1. — Reconstructed images of the disk seen face-on at z/=345GHz (left column) and 
900 GHz (right column). Upper row: Undisturbed disk, distance 50 pc. Middle/ lower row: 
Evolved disk with baroclinic instability, distance 50pc/140pc. The contour lines mark steps 
of 0.5mJy/beam and 1.5 mJy /beam at z/=345GHz and z/=900GHz, respectively. Total 
integration time: 2h. 
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Fig. 2. — Reconstructed images of the evolved disk seen under inclinations of i=0° (face- 
on; upper row), 30° (middle row), and 60° (lower row). Distances of 50 pc, 140 pc, and 
300 pc are considered. z/=900GHz; total integration time: 2h. The contour lines mark steps 
of 2.0 mJy/beam^ 1 '' and 1.5 mJy/beam 1 - 2 ^ in the case of d=50pc^ 1 ' ) , 140pc ( - 1 ^ and 300PC*- 2 -*, 
respectively. 
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Fig. 3. — Reconstructed images of the evolved disk seen face-on in distances of 50 pc and 
140 pc (z/=900GHz). Total integration times of 5min, 24min, 2h and 8h are assumed (from 
top to bottom). The contour lines mark steps of 2mJy/beam. 



